#/usr/bin/Rscript

on_100 <- read.table("on/profits.dat")
op_100 <- read.table("rf/profits.dat")
ut_100 <- read.table("ut/profits.dat")
iosup_100 <- read.table("iosup/profits.dat")
super_100 <- read.table("super/profits.dat")

print("100 op")
shapiro.test((op_100$V1 - on_100$V1) / abs(on_100$V1))
print("100 ut")
shapiro.test((ut_100$V1- on_100$V1) / abs(on_100$V1))
print("100 iosup")
shapiro.test((iosup_100$V1 - on_100$V1) / abs(on_100$V1))
print("100 super")
shapiro.test((super_100$V1- on_100$V1) / abs(on_100$V1))


#QQplots

jpeg("normalidade_g_100.jpg")
par(mfrow=c(2, 2))

qqnorm((op_100$V1 - on_100$V1) / abs(on_100$V1), main="op 100 0")
qqline((op_100$V1 - on_100$V1) / abs(on_100$V1))
qqnorm((ut_100$V1 - on_100$V1) / abs(on_100$V1), main="ut 100 0")
qqline((ut_100$V1 - on_100$V1) / abs(on_100$V1))
qqnorm((iosup_100$V1 - on_100$V1) / abs(on_100$V1), main="iosup 100 0")
qqline((iosup_100$V1 - on_100$V1) / abs(on_100$V1))
qqnorm((super_100$V1 - on_100$V1) / abs(on_100$V1), main="super 100 0")
qqline((super_100$V1 - on_100$V1) / abs(on_100$V1))


jpeg("homocedasticidade_g_100.jpg")
par(mfrow=c(2, 2))

plot((op_100$V1- on_100$V1) / abs(on_100$V1), main="op 100 0")
plot((ut_100$V1- on_100$V1) / abs(on_100$V1), main="ut 100 0")
plot((iosup_100$V1- on_100$V1) / abs(on_100$V1), main="op 100 0")
plot((super_100$V1- on_100$V1) / abs(on_100$V1), main="ut 100 0")

dev.off()



